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I. INTRODUCTION 



The parity operator V for a quantum field theory has the effect 

V(p{'x, t)V = X, t) (scalar field), 

P</p(x, t)V — — (/?(— X, t) (pseudoscalar field). 

We can also define an intrinsic parity operator Vj. This operator has the same effect as V 
except that it does not reverse the sign of the spatial argument of the quantum field: 

Vi(p{'x,t)Vi = (p{'x,t) (scalar field), 

Vi(p{'x., t)Vi — — <^(x, t) (pseudoscalar field). 

To determine the Lorentz transformation properties of an operator in a quantum field 
theory one must calculate the commutator of this operator with the generators of the Lorentz 
group Ji^^. If this operator does not depend on the space-time coordinates (x, i) and it 
commutes with J^" , then it is a scalar under Lorentz transformations. 

For the purposes of this paper we will suppose that the Hamiltonian H of our quantum 
field theory has parity symmetry; that is, [P, H] = 0. An example of such a theory is the 
free scalar quantum field theory whose Hamiltonian is 

j dx{|7r2(x,i) + |[V(/p(x,t)]V|/xV'(x,t)}. 
For such a theory it is easy to see that the intrinsic parity operator is a scalar because 

[Vi, J'^n = 0. 

However, while the conventional parity operator is a rotational scalar because it com- 
mutes with the three generators of spatial rotations, 

it does not commute with the Lorentz-boost generators: 

Therefore, the parity operator V is not a Lorentz scalar. Furthermore, as we will see, V is not 
the spin-0 component of a vector, a tensor, or indeed any finite-dimensional representation 
of the Lorentz group. 
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This observation raises a fundamental question: How does V transform under the Lorentz 
group? In this paper we provide a partial answer to this question. We argue that V 
transforms as an infinite direct sum of finite-dimensional tensorial representations of the 
Lorentz group. 

This paper is organized as follows: In Sec. Owe review the general theory of the irreducible 
representations of the Lorentz group. Then in Sec. IIIII we show that V transforms as a 
direct sum of irreducible representations of the Lorentz group, and that to identify each of 
these irreducible representations one must solve a difference-equation eigenvalue problem. 
In Sees. IIVI and E] we perform an analysis of special cases of this difference equation and 
show that the eigenfunctions are closely associated with the Wilson polynomials. Finally, 
in Sec. IVII we make some concluding remarks about the structure of the solutions to the 
general difference-equation eigenvalue problem that is derived in Sec. IIIII The properties of 
the Wilson polynomials are summarized in the Appendices. 



II. IRREDUCIBLE REPRESENTATIONS OF THE LORENTZ GROUP 

In this section we review briefly the classification of irreducible representations of the 
Lorentz group. Following the exposition of the irreducible representations of the Lorentz 
group by Gel'fand, Minlos, and Shapiro jlj. Bender and Griffiths and Griffiths we 
note that all irreducible representations of the Lorentz group are characterized and labeled 
by a pair of numbers (£o,^i)- The first number, Cq, is a nonnegative integer and the second 
number, £i, is in general complex. Any irreducible representation of the Lorentz group is 
a direct sum of irreducible representations of the rotation subgroup S'0(3) and £o is the 
lowest-spin of the representations contained in this sum. 

If £i — £o = or noninteger, then the representation belongs to the nonsingular class 
of infinite-dimensional representations. Such a representation contains an infinite tower of 
spins beginning with ^q! the spins £ in this representation are £ = + 1; + 2, ■ ■ ■ . 

Each spin in this sequence occurs once and only once. 

If — ^0 is a nonzero integer, then the representation belongs to the singular class. For 
this case if > io, the representation is finite-dimensional; it contains all spins i with 
^0 < ^ < l^il ~ 1 and each spin occurs exactly once. If < io, the representation is 
infinite-dimensional; it contains all spins i: io < i < oo and again each spin occurs exactly 
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once. 

Since the Lorentz group is noncompact, it also has indecomposable representations. How- 
ever, this situation does not arise in this paper. 

To describe the Lorentz transformation properties of a space-time-dependent quan- 
tity, such as a quantum field, we must specify all of the spin components of the field 
{Q^Qi Q^o+i' Qe.o+2i ■■■)• Let us suppose first that this field transforms irreducibly un- 
der the Lorentz group and let us also suppose that the lowest-spin component has spin 
0. Then, an infinitesimal Lorentz transformation of the spin-0 component Qo(x, t) has an 
orbital and a spin part. Specifically, 

-I [Qo, ^°'] = td'Q, - x'd'^Qo + Ql, (1) 

where Q\ is the spin-1 component of the irreducible representation. Under an infinitesimal 
Lorentz transformation, this spin-1 component then transforms like 

-I [Q{, J"*] = td'Qi - x'd'Qi + Q'i + «5*^go, (2) 

where is the spin-2 component of the representation. Note that is a pure spin-2 object 
and is therefore symmetric and traceless: Q2 — 0- (Repeated indices indicate summation.) 
The number a in (j^)) is directly related to the parameter ^l for the irreducible representation 
by the formula 

a = \{il-l). (3) 

This process of evaluating commutators can be used iteratively to calculate all of the spin 
components of the irreducible representation. 

Now suppose that Qo is the spin-0 component of a representation of the Lorentz group, 
but that this representation is not irreducible. This means that there may be more than one 
component of each spin. Thus, under an infinitesimal transformation we have 

-t [Qo, J°*] = td'Qo - x'd'Qo + iVl(Qo), (4) 

where is one of the spin-1 components of the representation. An infinitesimal Lorentz 
transformation of this spin-1 component produces a spin-2 component N2{Qo), for which 
^2iQo) = 0, and a new spin-0 component Nq{Qo): 

-I [NiiQo], = td'NiiQo) - x'd^NiiQo) + Nl^Qo) + 5''No{Qo). (5) 
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Since this representation is not irreducible, Nq^Qq) is not a multiple of Qq. 

This process of calculating commutators can again be used to generate all of the spin 
components of the representation. For example, commuting N2 {Qo) with the Lorentz boost 
will define a spin-3 component N^-'^{Qq). However, the procedure is significantly more com- 
plicated than that for an irreducible representation because there may be many components 
for each spin. Indeed, it may be that there are an infinite number of spin-0 components, Qo, 
No{Qo), Nq{Qo), and so on, and we will see that this is the case with the parity operator V. 

III. INFINITESIMAL LORENTZ TRANSFORMATION OF THE V OPERATOR 

Like the generators of space-time translations P'^, the parity operator P is a global 
symmetry operator and is not an explicit function of the space-time coordinates x^; that is, 
Qt^jy = 0. As a consequence, the orbital terms in (@)) vanish and we have 

Ni{V) = -i[P, J°'] = 2iJ^'V. (6) 

This identifies the spin-1 term that arises after commuting with the generator of a Lorentz 
boost. 

A second commutation with the Lorentz boost gives a spin-2 component and a new spin-0 
component: 

-z[Nl{V), J''] = Ni\V) + 5^^iVo(P), (7) 
where the new spin-0 component is 

iVo(P) = |J°V°'P. (8) 

Since Nq{V) is not a constant multiple of V, we conclude immediately that V does not 
belong to an irreducible representation of the Lorentz group. 

To analyze the Lorentz transformation properties of the parity operator it is convenient 
to introduce the following notation. First, we define two 3- vectors: 

V = \e'^^J^^. (9) 
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These operator quantities satisfy the following commutation relations: 



[K\K^] = -le'^'^L^. (10) 



Next, we define three operator products: 

W = K'K' 



A = UU 



m = K'V = VK\ (11) 

These three operators are mutually commuting: 

[W, A] = [A, m] = [m, W] = 0. (12) 

The W, A, and m operators have the following commutation relations with the generators 
of the Lorentz group: 

[W, K'] = [A, K'] = -2it'^^UK^ - 2K\ 

[W, U] = [A, U] = [m, K'] = [m, U] = 0. (13) 
Next, we define the operator B: 

B = A-W. (14) 
Observe that from (fT!?|l the operator B commutes with the generators of the Lorentz group: 

[B,K'] = [B,V] = 0. (15) 
Also, B commutes with the parity operator V: 

[B,V] = 0. (16) 

The vanishing of the commutation relations in (jl5|) and imply that B transforms as a 
scalar. The operator m also commutes with the generators of the Lorentz group but m is a 
pseudoscalar because 

VmV = -m. (17) 



Finally, we define the scalar operator M by 



M = m^. (18) 

The four operators B, M, W, and parity V are a mutually commuting set: 

[B, M] = [B, W] = [M, W] = [B, V] = [M, V] = [W, V] = 0. (19) 

We know that the parity operator does not transform irreducibly under the Lorentz group. 
Thus, we will assume that it transforms as a direct sum of operators that do transform 
irreducibly. From the commutation properties of the operators defined above we know that 
the most general operator that could be produced by repeated commutation of V with the 
generators of the Lorentz group has the form f{B, M, W)V, where f{B, M, W) is an as yet 
unknown function of the operators B, M, and W. We can see from (jHJ that this is the most 
general spin-0 structure that can appear with repeated commutation. Let us assume that 
the function f{B, M, W) has a power series representation of the form 

f{B, M,W)=Y, ai^nB'M^'W''. (20) 

l,m,,n 

A lengthy calculation of the double commutator of f{B, M, W)V with the generators of 
the Lorentz group gives the new spin-0 component 

No[f{B, M, W)V] = -l(w + ^j^] f{B, M, W)V ^ 



3\ B + WJ'' ' 3V1 + 4:B + 4W 

2B + 4W+(w- ^^y^ ^ (Vl + 45 + 4iy - l) 

X f (^B,M,W + 1 + Vl + 45 + AW^ 

-2B-4W+ (w- (Vl + iB + 4W + l 

V B + W J \ 

xf (b,M,W + 1- Vl + ^B + 4W^^V. 



(21) 



Let us now suppose that f{B,M,W)V transforms as the irreducible (0,£i) representation 
of Lorentz group. Then, from (jH)) we have 

iVo[/(5, M, W)V] = Uil- l)/(5, M, W)V. (22) 
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Combining (^1]) and (j^ . we conclude that f{B,M,W) satisfies the functional relation 



This functional equation is a difference- equation eigenvalue problem in which the eigen- 
value (1 — if) determines the representation of the Lorentz group under which the operator 
f{B, M, W)V transforms. We do not yet know how to find the general analytical solution 
to this equation. However, in the next two sections we show how to solve analytically this 
master equation for two important special cases. 

IV. SPECIAL CASE I: S = AND M = 

In this section we solve an important special case of the difference-equation eigenvalue 
problem (j2Sl)- Let us assume that 5 = and M = 0. Then, pHjl simplifies to 




x/ [B,M,W + 1- Vl + ^B + AW 



')\ = {l-il)f{B,M,W). (23) 




(3 - v^TTW) W 



f (w +1- Vl + 4W^) = (1 - iDfiW). (24) 



To solve this equation we make the change of independent variable 




(25) 





i). Next, we introduce the dependent variable 



(26) 



The functional equation satisfied by g{z) is 



{2z + l){2z + 3fg{z + 1) - 4^(2^ - l){2z + l)g{z) + {2z - l){2z - 3fg(z - 1) 
8z{il-l)g{z). 



(27) 
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In order to determine the eigenvalues (il — l) we must impose an eigenvalue condition. The 
condition that we will use is that f{W) be an entire function of W. This choice is a natural 
physical constraint because it eliminates the possibility that the irreducible representation 
of the Lorentz group in (j22j) could be a nonsmooth function of the physical parameters B, 
M, and W. A similar eigenvalue condition based on the entirety of a solution to a difference 

n 

equation may be found in Ref. j^. 

Imposing this constraint gives an infinite sequence of allowed values for ii: 

£i{n)=2n + l (n = 0, 1, 2, 3, ■ ■ ■ )• (28) 

Note that the eigenvalues involve the square of ii but we may assume that ii is positive . 
The irreducible representations of the Lorentz group corresponding to these values of ii{n) 
are traceless, totally symmetric tensors of rank 2n: T, T^'^ , rp/iuXa ^ Thus, we 

see that the parity operator transforms as a direct sum of the spin-0 components of these 
tensors [see (HUj)]. 

The eigenfunctions corresponding to the above eigenvalues are all polynomials in 2;^ except 
for the eigenfunction corresponding to the lowest eigenvalue for which the eigenfunction is 
a rational function: 

1 



9oiz) 



1 1 ' 
7 — — 

^ 4 



g^{z) = + |, 

9^{z) = / + \z' + If, 

g,{z) = + fz' + + 

g,{z) = / + 19/ + '-fz' + '-^z' + (29) 

and so on. We have normalized the eigenfunctions so that each of the polynomials is moniq 
that is, the coefficient of the highest power of z"^ is unity. 

The polynomial eigenfunctions in (j^^ are Wilson polynomials W„(x^; a, 6, c, ci), where 
the parameters are given bya = 6=l/2 and c = d = 3/2 (see Appendix A): 

^^n(.) = ^^W„_i(-.^;l,i,|,f). (30) 

These polynomials can be expressed as generalized hypergeometric functions: 

9n{z) = ^ ^-j^ ^4F3 (1 - n, n + 2, i - ^, i + z; 1, 2, 2; l) . (31) 



From the solutions go{z) in and gn{z) in (|HT|) we can construct the eigenfunction 
solutions to 



fn{W) 



i-V^+y^ (n = 0), 



(n--l)!(n!)2(n+l)! (^32) 



X 4F3 (l - n, n + 2, i - + i, i + + i; 1, 2, 2; l) (n > 0). 
Here are the first six solutions: 



/^(ly) = e^-V^+V4vf/ ^ 20Vr3 ^ ^Qgyf^2 ^ ^7gyj^ ^ 4|o^ _ (33) 

The difference equation ()27p is second order and linear, and this means that for each 
eigenvalue there are two linearly independent solutions. It is straightforward to use the 
method of reduction of order (which is ordinarily used for linear differential equations) to 
calculate the second solution. We seek a second solution of the general form 

hn{z) = gn{z)Un{z), (34) 

where m„(z) is an unknown function to be determined. The function Un{z) is easy to find 
because it satisfies a first-order difference equation: 

Un{z + 1) - Un{z) {2z - l){2z - ?>fgn{z - 1) 



Un{z) - Un{z - 1) (2z + l)(2z + 3) + 1) ' 

From this equation we find that apart from a multiplicative constant, we have 



(35) 



Un{z) - Un{z - 1) = — ^.3.^^ , , 7^^- (36) 

{2z - 3)^(22; - iy\2z + lygniz - l)gn{z) 
Summing both sides of this equation, we get Un{z) (apart from an additive constant): 

^ 1 
Un{z) = (2a; - 3)2(2x - l f{2x + 1) - l)gn{x) ' ^^^^ 

X — Zq 
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where zq is arbitrary. 

This sum can be evaluated analytically for the case n = and we obtain 

hoiz) = ^-^1^. (38) 

When n > we must leave the second solution in the form of a sum: 

^ 1 
h^iz) = gn{z) J2 (2x - 3)2(2x - l)3(2x + l^g^ix - l)g4x) " ^^^^ 



x=zo 



The general solution to ()27|1 is a linear combination of gn{z) and hn{z). However, as a 
quantization condition, if we demand that the solution to ()24|) be entire, then the coefficient 
of hn{z) must vanish. 

Our objective now is to reconstruct the parity operator as a direct sum over the irreducible 
representations constructed by multiplying fn(W) in (jH^ by the parity operator V: 

oo 

V = J2<^nfn{W)V, (40) 

n=0 

where {c„} are coefficients that must be determined from the equation 



Co + (^^ - i) X] ^"S-nl^) 

n=l 



(41) 



l = 5Zc„/„(iy)=e^ 

ra=0 

which is obtained from ()40j) by multiplying by V. 

Equation (PT|) must hold for every value of z. Thus, if we let z = |, we obtain 

Co = -t. (42) 

For n > 1, the coefficients c„ satisfy 



z^ 



T- = Y.^ng^n\^) = $^(-l)"c.+iP.(-Z^), (43) 

4 n=l n=0 

where P*„(x^) is defined in (IA1|) . We solve this formal functional equation by continuing 
it analytically into the complex 2;-plane and then making the change of variable x = iz. 
Equation (jlHjl then becomes 



4 (e"^'^ + i) 



1 + 4x2 

n=0 



5^(-l)"c„+iP„(x2). (44) 
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We then multiply both sides of this equation by j7r'^x{l + Ax'^y smh.{x7r) Pm{x'^) / cosh^(x7r) 
and integrate over x from to oo. Using the orthogonality property of the Wilson polyno- 
mials in ()A2|1 . we obtain the following quadrature formula for c„ for n > 1: 

2 (-l)"[(2n)!]2(2n + l)! . (e-^- + z) sinh(x7r) ^ . ^. 

c- = ^ Tr iMi2/n4r^ , 1M14 / dxx{l + Ax^y ' ' Pn-iix"). 45 

[(n - + 1)!]4 Jo cosh'^(x7r) 

Note that this analytic continuation is an extremely delicate process because the series 
in ()40|) has only a formal existence. Indeed, while the Wilson polynomials ()A1|) are com- 
plete and orthogonal, the polynomials obtained by replacing by —z^ [see fHH]) ] are not 
orthogonal with a positive weight function and not complete in the usual sense. Thus, the 
analytic continuation above is a procedure that converts a purely formal series identity into 
a series that actually converges with coefficients that can be determined. 

If we substitute the expression for c„ {n > 1) in ()45|1 and also the value of cq in ()42|l 
into ()40|). then ()40|) becomes an identity. This implies that we have decomposed the rep- 
resentation under which V transforms into a direct sum of finite-dimensional irreducible 
representations 

(0, 1) © (0, 3) © (0, 5) © (0, 7) © ■ ■ ■ . (46) 

That is, we have shown that V transforms as a scalar plus the spin-0 component of a 
two-index symmetric traceless tensor plus the spin-0 component of a four-index symmetric 
traceless tensor plus the spin-0 component of a six-index symmetric traceless tensor, and so 
on. This is the central result of our analysis. 

We conclude this section with the remark that the first term in the direct sum above is a 
scalar; that is, the term fo{W)V transforms like a Lorentz scalar. Thus, V transforms like 
the elementary structure q-^'^V^^^/^. 

V. SPECIAL CASE II: M = 

In this section we solve an eigenvalue equation that is much more general than that solved 
in Sec. IIVI We consider here the case M = but we allow B and W to be arbitrary. Now, 
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simplifies to 

2Wf{B,W) + 



X 



(2B + 3W + WVl + 45 + 4W^ f (b,W + 1 + Vl + 45 + 4W 



2B + 3W- WVl + AB + AW) f (B,W + 1 - + 4:B + AW 



(47) 



To analyze this equation we generalize the substitution that we made in (|^^ and make 
the change of independent variable 



W + B 



(48) 



In terms of this new variable we have f{B, W) = f (5, 2;^ — 5 — |) . Next, we introduce the 
new dependent variable by the substitution 



(7(5,^)=e— 7 {B,z'-B~\) 



(49) 



which is a generalization of (|^. We obtain the following functional equation satisfied by 

g{B,zy. 



2{z'-B--)g{B,z)-- 



3z' - B - - + 2z { - B - - 



4 



9iB,z + l) 



1 

+2I 



3z' -B 2z U - 5 - - 



4 



4 



g{B,z-l) = {l~£l)g{B,z). 



(50) 



The solutions to this difference equation are Wilson polynomials W„(x^; a, 6, c, c/) with 
parameters a = 6 = |,c=| — ^/B + 1, = | + ^/ B + 1 (see Appendix B). Thus, 



9n{B,z) 



{2ny. 



-z^, i |, i - VbTT, \ + v/5TT 



{n\fT (n + 1 - VbTI) T (n + 1 + /bTT) 



(2n)! r (1 - VeTI) r (1 + VSTT) 

X 4F3 (-n, n + 1, 1 - 2, i + 2; 1, 1 - + 1, 1 + V5 + 1; 1) . (51) 
The first four of these polynomials are 



go{B,z) = 1, 

g,{B,z) = z'-lB-l 

g2{B,z) = z'-{B-l)z' + lB'-lB-^^, 

g,iB, z) = z^- (IB - f ) + ilB^ -KB + g)z^- ±B^ + - f^S - ill. (52) 
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Correspondingly, the solutions of (jTTjl are 

^'^^ (2r^)!^(l-^/En)^(l + V5TT) 

X 4F3 (^n, n + 1, i - + i + ^Jw+B + l; 1, 1 - ^/B + l, 1 + V^+T; 1^ (53) 



and the first four of these solutions are 



/3(i?, ly) = e^V^+^+V4 ^^^3 + 52 (Spy + f ) + 5 + f ly + I) 

+iy (iy2 + 4iy + . (54) 

It is interesting that while we are now treating B as an arbitrary parameter, the eigenvalues 
are independent of B and thus are the same as in the case considered in Sec. |TV| [see (j2HI)]- 
To find the Lorentz transformation properties of the parity operator we need to recon- 
struct the parity operator in terms of these solutions fn{B, W): 

00 

V = Y,Cnfn{B,W)V. (55) 

n=0 

This requires that we find the coefficients c„ in the identity 

00 00 
1 = ^«^-(^' ^) = ""-SniB, z). (56) 

n=0 n=0 

To analyze this identity we multiply it by e~*'^^ and obtain 

00 

e-^-^ = J2cn9n{B,z). (57) 

n=0 

Next, we continue analytically this formal functional equation to the complex z-plane and 
let X = iz. Equation (|57j) now becomes 

00 

e-"" = 5^(-l)"c,P„(a;'), (58) 

n=0 

where Pnix"^) is defined in ()B1|) . Next, we multiply the both sides of (fHS|) by 

47r2(2n)! {2n + 1)! x tanh(a;7r) 



(n!)4r2 (n + 1 - y/B~+T) {n + 1 + ^/BTT} (cos {2n^/B~+T) + cosh(2x7r)) 
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and then integrate over x from to oo. Using the orthogonahty property of Wilson polyno- 
mials in ()B2|) . we get the formula for Cn- 



We conclude that we have decomposed the representation under which V transforms into 
a direct sum of finite-dimensional irreducible representations. Even though the operator B 
is taken to be nonzero, the eigenvalues in the difference-equation eigenvalue problem ()47|) 
remain unchanged. Thus, the conclusion of Sec. IIVI that the parity operator V transforms 
as the direct sum of the spin-0 components of the finite-dimensional tensor representations 
(0, 1) © (0, 3) © (0, 5) © (0, 7) © ■ ■ ■ remains unchanged. 

Finally, as we observed at the end of Sec. IIVI we point out that the first irreducible rep- 
resentation /o(-B, W)V in the direct sum is a scalar. Thus, we may conclude that under 
a Lorentz transformation the parity operator V transforms as the operator g-*'rV'S+M/+i/4^ 

VI. FINAL REMARKS 

The general case for ()23|) is obtained when all three parameters, B, M, and W, are 
nonzero. For this case we cannot solve ()23]) in terms of polynomials but we conjecture that 
the eigenvalues in ()28|) remain unchanged [6|. Assuming that this is indeed the case, we may 
conclude that the Lorentz transformation properties of the parity operator are unchanged 
from what we found in Sees. IIVI and IV j namely, the parity operator transforms like the direct 
sum of irreducible representations in (jlUj) . 
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4(-l)'^7r2(2n)! (2n+ 1)! 



-n 



{niyr^ {n + i- VbTT) r2 (n + 1 + VbTT) 




(59) 
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APPENDIX A: PROPERTIES OF WILSON POLYNOMIALS 



In this appendix, we list the properties of the Wilson polynomials W„(x^; a, b, c, d) with 
parameters a = 6 = |,c = (i=|0,0. For simplicity, we define the associated monic 
Wilson polynomials 



IW (x'^-l 1 3 3\ 
^'^"V '2'2'2'2/' 



(Al) 



(2n + 2)! 

(The coefficient of the highest power in a Monic polynomial is unity.) These polynomials 
are used in Sec. IIVI 

The orthogonality condition is 

^2 sinMxvr) „ , 2^ o . 2^ (n!)2[(n + l)!]^[(n + 2)!]^ 



vr 



— / dxx(l + Ax 



2\2 



cosh (X7l) 



[(2n + 2)!]2(2n + 3)! 



-5n 



The polynomials P„(a;2) satisfy a three-term recurrence relation: 



PJx' 



x^ — -n{n + 1) 



Pn-lix^)+ 



where the first two polynomials are Pq{x^) = 1 and Pi{x'^ 
Some generating functions for these polynomials are 

2(2n + 2)!P„(a;2)t'^ 



in - lYn\n + 1)^ 
4(2n-l)(2n + l) 

2 3 



■Pn-2{X') 



X 



4" 



n=0 
oo 



(n!)2[(n + 2)!]2 

(2n + 2)!P„(a;2)r 
^n![(n + l)!p(n + 2)! 

(2n + 2)!P„(x2)r 



E 



^Fi (i + zx, i + zx; 1; -t) 2F1 (| - ia;, | - zx; 1; -t) 
iFi (I + zx, I + zx; 2; -t) 2F1 (| - ix, | - zx; 2; -t) 



E 

n=0 



(n!)2[(n + l)!]^ 



r4i'3 



4t 



2, 2, 2 + ix, 2 ^a;; 1, 2, 2; ^ ^^2 



(A2) 



(A3) 



(A4) 



APPENDIX B: FURTHER PROPERTIES OF WILSON POLYNOMIALS 



In this appendix, we list the properties of the Wilson polynomials W„(x2; a, 6, c, d) with 

1 



parameters a = b - 
monic polynomials 



2' 2 



v/pTI, = i + VpTT 



We define the associated 



PJx' 



{2n)\ H '2' 2' 2 



VP+T,i + VP+T . 



(Bl) 



These polynomials are used in Sec. |3 
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Orthogonality: 



47r' 



tanh(x7r) 



cos (27r/BTT) + cosh(2 



XTT] 



[n 



(2n)! (2n + 1)! 



Recurrence relation: 

, B 1 + 



PJx' 



^2 4: 



E 



^ (n!)2r (n + l - /bTT) r (n + 1 + VeTT) 

sin 



, ^2Fi {\ + ix,\- V5 + 1 + zx; 1 - + 1; -t) 

X 2F1 -ix,\ \ VB + 1 - ia;; 1 + V5 + 1; -t 

sin 

7r(l + t)v/5^ 



=4F3 (^|, 1, i + ^a;, i - tx; 1, 1 - V5 + 1, 1 + VB + 1; (^^^^ 



(B2) 



(B3) 



4(2n- l)(2n + 1) 

with Po{x^) = 1 and Pi{x^) = + 1^ + 1. 

Generating functions: 

^ (2n)!P„(x2)r 
^ (n!)4 

n=0 ^ ^ 

= 2F1 (i + ix,l + ix; 1; -t) 2F1 (^^ - + 1 - ix, i + + 1 - tx; 1; -t^ (B4) 

and 

(2n)!P„(x2)r 



(B5) 
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